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£NJ ' Abstract 

^— ^ ' We consider the energy-efficiency of Multiple-Input Multiple-Output (MIMO) systems with constrained received power rather 
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than constrained transmit power. A Energy-Efficient Water- Filling (EEWF) algorithm that maximizes the ratio of the transmission 
rate to the total transmit power has been derived. The EEWF power allocation policy establishes a trade-off between the transmission 
rate and the total transmit power under the total receive power constraint. The static and the uncorrelated fast fading Rayleigh 
channels have been considered, where the maximization is performed on the instantaneous realization of the channel assuming 
perfect information at both the transmitter and the receiver with equal number of antennas. We show, based on Monte Carlo 
simulations that the energy-efficiency provided by the EEWF algorithm can be more than an order of magnitude greater than 
the energy-efficiency corresponding to capacity achieving Water-Filling (WF) algorithm. We also show that the energy-efficiency 
increases with both the number of antennas and the signal-to-noise ratio. The corresponding transmission rate also increases but 
at a slower rate than the Shannon capacity, while the corresponding total transmit power decreases with the number of antennas. 

c/j ' I. Introduction 

Energy consumption is a key issue in the deployment of green wireless communication networks [1]. Currently, there is a 
marked increase in the number of multimedia functionalities available at the mobile terminal that require ever higher transmission 
^ rates. This, together with the large screen size of mobile terminals, directly translates into a higher power consumption. On 
, the network side, the energy-efficiency depends, among other things, on the power transmitted at the base station. Higher 
£SJ ■ transmission rates require more base stations to cope with the increase in data traffic. However, the gains resulting from heavy 
0^ , cell splitting tend to be severely limited by high inter-cell interference. Moreover, high CAPEX (capital expenditures) as well as 
■ high OPEX (operating expenditures) associated with high power macro nodes further limits the usefulness of such an approach 

in : [2]. 

Within this context, the introduction of the Multiple-Input Multiple-Output (MIMO) concept has become instrumental to 
achieve high spectral and energy-efficiencies [3]. MIMO employs multiple antenna elements at both the transmitter and the 
receiver and has attracted extensive attention due to the potential manifold increase in performance it can bring to standardized 
wireless communication networks, such as, the UMTS/LTE/LTE-A, WIMAX and WIFI systems as well as systems beyond 
these. Indeed, the multiple antennas can be exploited by creating a highly effective antenna diversity system which combats 
$J] \ the effects of fading so as to improve the average signal-to-noise ratio (SNR) or, by spatial multiplexing of the transmit data, 
& ■ i.e., the data streams are spatially separated by the multiple antennas at the receiver and the transmitter [4]. 

The energy-efficiency analysis of MIMO systems has attracted the attention of the communication research community 
lately. A survey on energy-efficient communications can be found in [5]. As outlined there, two main research approaches can 
be identified:(l) the pragmatic approach, which focuses on system specific features such as modulation and coding-decoding 
schemes as well as electronics and, (2) the information-theoretic approach, which focuses on the maximization of channel 
capacity per unit cost function. The latter approach has been considered in [6] where the energy-efficiency is defined as the 
ratio between the transmission rate and the transmit power. They show that the optimum energy-efficiency is obtained for 
transmit power that tends to zero for static and fast-varying MIMO channels with informed transmitter and receiver. 1 This 
result illustrates the fact that according to the provided energy-efficiency definition the maximum is achieved by transmitting at 
very low powers. At the same time this requires the system to support very low data rates as, for example, in sensor-networks 
[7-9]. However, as they noted in [6], this is not a practical result for communication systems where a given data rate is required 
for satisfactory operation of the wireless network. Hence, a different strategy is required to obtain a non-trivial solution to the 
energy-efficient communication over MIMO wireless channels. This is what we investigate in this article. 

The solution we have found is based on the following observations. Firstly, limiting the maximum total transmit power 
while maximizing the energy-efficiency, as we have seen, is not worthwhile. Secondly, we find that the transmit power can 
still be constrained, though indirectly, by imposing a constraint on the total receive power. This constraint is justified by the 
fact, that there is a minimum power level (or signal-to-noise ratio) required at the receiver for a satisfactory operation of the 
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'For the slow-fading MIMO channel they show that a non-trivial solution exists, i.e., non-zero transmit power allocation scheme. 
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communication link. Below this minimum level, which corresponds to the receiver sensitivity, the transmit signals cannot be 
decoded properly which leads to unwanted performance degradation. It turns out that this minimum receive power constraint 
can be used to obtain the power allocation for a MIMO channel that maximizes energy-efficiency. 

In this paper we adhere to the information-theoretic approach outlined in [6] and find the optimal strategy for energy-efficient 
transition satisfying a total receive power requirement. We specialize our analysis to the static and the fast-varying single-user 
MIMO channels leaving the slow varying MIMO channel to be considered elsewhere. We also assume that channel state 
information is available at both the transmitter and the receiver. The contributions of this paper are summarized below : 

• We derive a Energy-Efficient Water-Filling (EEWF) algorithm that maximizes the ratio of the transmission rate to the 
total transmit power under total receive power constraint. This is a novel approach and differs from the customary total 
transmit power constraint that gives the trivial zero power result [6]. 

• We show that the static isotropic 2 and the static rank-1 MIMO channels achieve the Shannon capacity if they are power 
efficient in the sense of the EEWF power allocation algorithm. This is obtained under the assumption of equality of 
the signal-to-noise ratios used in the EEWF and Water-Filling (WF). Here we show that while the optimum transmit 
power goes to zero as the number of antennas increases, the transmission rate also increases but to a limit defined by the 
minimum total received power and the noise variance. 

> We show that static SISO channels also achieve the Shannon capacity under the conditions given above; however, for the 
Rayleigh fading SISO channel this is only observed in the high signal-to-noise ratio limit. 

• We derive an upper and a lower bound for the ratio between energy-efficiency corresponding to the EEWF and the WF 
power allocations as well as the ratio of the transmission rate corresponding to the EEWF and the Shannon SISO capacity 
under the assumption of finite inverse transmit power. 

• We illustrate for Rayleigh fading MIMO channels, based on Monte Carlo simulations, that the energy-efficiency provided 
by the EEWF algorithm can be more than an order of magnitude greater than the energy-efficiency corresponding to 
capacity achieving Water-Filling (WF) algorithm; while the corresponding total transmit power decreases with the number 
of antennas. 



II. Static MIMO Channel 

The input-output relationship for a MIMO communication system with N t transmit antennas and iV r receive antennas can 
be written as 

y = Hx + n, (1) 

where y 6 C^xi is the received signal vector, x 6 C tXl is the transmit signal vector, H 6 C tX * is the MIMO channel 
matrix and n 6 C NrXl is the noise (AWGN) vector with covariance matrix lZ n = (nn*) = Ogl/v,., where ctq is the noise 
variance, (} denotes expectation and (y denotes the Hermitian transpose operation. 
We define the total transmit power as 

if = ti{n x }, (2) 

where 1Z X = (xx*) is the covariance matrix of the transmit vector and tr(A) denotes trace operation which is the sum of the 
diagonal elements of matrix A. 

The total receive power* is obtained by fixing the channel matrix H and is defined for the noiseless case as 

If = tr{HtH7^}, (3) 

As shown in [3], (1) can be transformed into a stream of r = min{Af r , Nt} parallel channels as follows 

y = A 1/2 x + n, (4) 

where y = U+y, x = V+x, and n = U+n. The matrices U E C N * xN * and V e C NtXNt are unitary and A 1 / 2 = diag{A,} /2 } g 
M. N ' xNt is non-negative and diagonal containing the singular values of the MIMO channel matrix, i.e. 

H = UA 1/2 V t . (5) 

Hence, the transmission rate of the MIMO channel (4) equals the transmission rate of the original channel (1) due to the 
unitary affine transformation and is given by 

R x = log 2 det(Ir + cr^A 1 / 2 ?^ 1 / 2 ), (6) 

We measure the energy-efficiency as the ratio between the transmission rate and the transmit power (bps/Hz/Watt) [6]. Hence, 
we divide (6) by (2) and obtain 

x _ Rf log 2 dct(I r + ( 7- 2 A 1 / 2 ^A 1 / 2 ) 

11 if tx{n x } ' Kn 

2 An isotropic channel is defined as a channel for which the likelihood of receiving a signal is the same for all directions of the unit sphere. 

3 It is worthwhile to mention here that the pathloss effects can be accounted for by multiplying the channel matrix with the corresponding factor. 
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We now formulate our energy-efficiency maximization problem with constrained receive power as follows 

log 2 det(I r + cr,7 2 A 1 / 2 72. x A 1 /2) 



(8) 



tr{n x } 

s.t. trjAV^A 1 / 2 } > P r , (9) 

where we have used (3) with P® > P r to get the constraint, P r is the minimum power required at the receiver and we use 

the MIMO channel power normalization tr{A} = N r N t . 

We further specialize our analysis to the case of equal number of transmit and receive antennas, i.e., N = N T = iVf 
Recalling that the rate (6) is maximized for independent Gaussian zero-mean complex transmit signal vector x with diagonal 

transmit covariance matrix [3] 

Tl x = diagfe}, (10) 

where pi > is the power allocated to channel i. 

We can now recast (8) and (9) into the following form 

max ELlog 2 (l + ^ 2 A^) ; (n) 
s.t. £*=i XiPi = Pri (12) 

Pi > o, 

where the power channel normalization is J2iLi = N 2 . 

The optimal transmit powers p = (pi,p2,--- ,Pn) can now be obtained by solving a nonlinear programming (NLP) 
optimization problem defined by (11) and (12). The solution must satisfy the Karush-Kuhn-Tucker (KKT) conditions [10]. We 
first construct the cost function containing the Lagrange multipliers v = (yi, l>2, ■ ■ ■ ,vn) f° r the inequality constraints and a 
multiplier /i for the equality condition 

^CP,",M) = (13) 



The KKT conditions that maximize (13) are 




V p J-(p,iy, M ) = ) (14) 
p}j = 0, (15) 

(16) 
(17) 
(18) 




where V p = (^-, gfj, • • • , gf^-) is the gradient operator. 
Thus, we get from "(14) that 



^5fef -Ef =1 ln(l + . - 2 « 



uXi + Vi = 0. (19) 
ln(2)(Et 1 P 4 ) 2 

Noticing that the constraint pj > need not be tight we can eliminate the constants since they act as slack variables. Hence, 
from (19) and (18) we get 

< V x - pif Ai, (20) 



ln(2)(^+ Pi ) 

where rj x , P t x , Aj and ctq have been defined above i = 1, , N. If — ^Pt^i < i^lfo 3 t l ien (^-* if K > 0- Hence, we 

get from (17), (18) and (19) that p q = ln ( 2 )( v -:-f J .p^\ i ) ~ xf • Now ' if 7 f ~ ^ p t^i > \ n (2)^ then Pi = since Pi > cannot 
satisfy (20) and (17) at the same time. We further see that 5Z i=1 \Pi = P r , i-e-, this constraint is tight which leads to fi > 
from (15). We have proved the following Proposition. 

Proposition 1 (Energy-Efficient Water-Filling (EEWF)): Consider the energy-efficiency of MIMO communication systems 
operating over a static AWGN channel and constrained minimum total received power defined according to (8) and (9). Then, 
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for informed transmitter and receiver the optimum is achieved by a "water-filling" power allocation solution pi with per-channel 
variable power levels according to 

Pr = ^(m^-^)- 02 ) ' (22) 

V X = ^F. (23) 
fl* = ^log^l + ^Aift), (24) 

i=l 

ff = (25) 

i=l 

where /i > is chosen to meet the receive power constraint (22), rf is the maximum energy-efficiency and R x and are the 
corresponding transmission rate and transmit power, respectively. The summation is over N x < N channels since the number 
of actual channels satisfying (21) and (22) may be less than or equal to N, Here (a) + denotes max(0,a) and ln() denotes 
the natural logarithm. 

As we can see, the optimization problem at hand involves finding p$, rf and /i. Finding n from (22) leads to finding the 
roots of a polynomial of at most degree N x in /i. Moreover, the solution is given by the largest positive real root as we show 
in Appendix A. 

It is worthwhile to note that the transmission rate corresponding to the maximum energy-efficiency always satisfies R x < C, 
where C is the channel capacity of the static channel achieved by the water-filling (WF) algorithm, [3]. In order to be able 
to compare different power allocation strategies we need to define a common reference. In our case we use the transmit 
signal-to-noise power SNR t = ~^P- . 

Corollary 1.1: The static SISO channel achieves the Shannon capacity if it is energy-efficient in the sense of Proposition 
1 and it satisfies the condition of equal transmit signal-to-noise ratio such that = P t , where, Ai is the channel gain. 

The above result is straightforward. Indeed, set A = 1 into (21)-(25) and compare it with the capacity of the AWGN SISO 

channel [3] for the same transmit signal-to-noise ratio SNR t = —5-. We see then that the EEWF power allocation policy(clearly 

a o 

there is only one channel) gives p\ = P r /\\. Hence, Pf = P r /Ai, if = Ailog 2 (l + Sjt)/P r and R x = log 2 (l + ^i). 

CT CT 

On the other hand, the SISO channel capacity is given by C = log 2 (l + -^j 1 ) with p\ = P t and P x = Pt\\- Hence, 
r] x c = log 2 (l + ^ L )/P t - Hence, R x = C when ^ = P t . 

This result has the following interpretation (static SISO case): If the transmit and receive power constraints satisfy = P t , 
then an increase or decrease of the transmit power will proportionally affect the optimum power allocation into the single 
link. This is not true for the general MIMO case, since the optimal power allocation in the sense of Proposition 1 and that 
achieving the MIMO capacity are not equivalent. Next we exemplify this latter statement by specializing our results to two 
special cases. We also look at their asymptotic behavior. 

Example 1: The Isotropic Static MIMO Channel. 

Under the isotropy assumption of the static channel all the eigenvalues of the MIMO channel matrix are equal. Hence, the 
normalization J2l=i ^« = N 2 gives Ai = N for i = 1,...,N, By further substitution into (21)-(25), we readily obtain the 
optimal transmit power allocation per channel that achieve maximum energy-efficiency and the corresponding total transmit 
power and transmission rate: 

Pi - ^, (26) 
N 2 P 

t = ^log 2 (l + -^), (27) 

R x = jVlog 2 (l + -p-), (28) 
erg A 

Pt = (29) 

Now compare the above with the capacity achieving strategy (WF) for the isotropic static MIMO channel with transmit power 
constraint ^=1 Pi = p ^ PL We have then that p t = rf = Alog 2 (l + ^)/P t , R x = C = Alog 2 (l + P x = NP f . 
Hence, we see that energy-efficient channel also achieves the Shannon capacity under certain conditions as stated in the 
following result. 



5 



Corollary 1.2: The isotropic static MIMO channel achieves the Shannon capacity if it is energy-efficient in the sense of 
Proposition 1 and it satisfies the conditions of equal transmit signal-to-noise ratio and equal number of antennas such that 

AT - *t- 

Remark 1.2.1: A energy-efficient, in the sense of Proposition 1, isotropic static MIMO channel with a large number of 
antennas behaves as follows 

lim Pi = 0, (30) 

7V-S- + 0O 

lim rf = oo, (31) 

N—^+oc 



_P 

2V->+ao" ln(2)er ( 



lim R x = W J 2 , (32) 



lim Ff = 0, (33) 

As we can see, the energy-efficiency increases asymptotically to infinity while the corresponding transmit power goes to zero 
and the transmit rate goes to a maximum limiting value. The physical interpretation is that if we constrain the receive power 
to a certain level, then the energy-efficiency will increase as we increase the number of antennas while the needed transmit 
power will decrease accordingly. 
Example 2: The Rank-1 Static MIMO Channel. 

Under the rank-1 assumption of the static channel, all the eigenvalues of the MIMO channel matrix but one are equal zero. 
Hence, the normalization J2iLi ^« = ^ 2 gi yes M = N 2 and = for i = 2, . . . , N. Further substitution into (21)-(25) gives 
the following results 

Pi = ^, (34) 
N 2 P 

rf = p- log 2 (l + -£), (35) 
R x = log 2 (l + % (36) 



T = (37) 

Corollary 1.3: The rank-1 static MIMO channel achieves the Shannon capacity if it is energy-efficient in the sense of 
Proposition 1 and it satisfies the conditions of equal transmit signal-to-noise ratio and equal number of antennas such that 

This result follows by noticing that the capacity achieving strategy for the rank-1 static MIMO channel with transmit power 
constraint ^2 i=1 Pi = Pt [3], is equivalent to the transmission rate and transmit power that corresponds to the maximum 
energy-efficiency and satisfying conditions above. Indeed, we have for the capacity achieving strategy that Pi = and 

~t- 

Remark 1.3.1: A energy-efficient, in the sense of Proposition 1, rank-1 static MIMO channel with a large number of 
antennas behaves as follows 



Pi = for i = 2, . . . , N, rf = log 2 (l + P t , R x = C = log 2 (l + P? = N 2 P t . 



lim p 1 = 0, (38) 

N-S.+00 

lim if = oo, (39) 

W-5- + 00 

lim R x = log a (l + -I), (40) 

lim Pf = 0, (41) 

As we can see, the energy-efficiency increases asymptotically to infinity at the same time as the corresponding transmit power 

goes to zero and the transmit rate remains constant and equal to the capacity of the SISO channel and is independent of the 
number of antennas. 



III. Fast Varying MIMO Channel 

We now consider a channel that undergoes rapid fading variations and is ergodic under the transmission duration. We further 
assume that the elements of the channel matrix H in (1) are now i.i.d. zero-mean complex Gaussian variables, i.e., the channel 
is Rayleigh fading. 

We analyze the energy-efficiency of a time-varying ergodic MIMO channel as the maximization of the instantaneous energy- 
efficiency given by (7) with an instantaneous constraint on the total receive power. The latter basically defines a power control 
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constraint on the link. Hence, we define the corresponding optimization problem as follows 

max^ 1062 ^^), (42) 

Ei=l Pi ' 

s.t. ^1=1 } iVi = Pr ' (43) 
Pi > 0, 

where the expectation () is taken over the realization of Ai and the maximization is over p±. The power channel normalization 

We can now construct a cost function and KKT conditions similar to (13) and (14)-(18), respectively. We enforce this 
conditions on the instantaneous realizations of A; and p t . Clearly, we seek find the instantaneous pi that maximizes the 
function within the expectation brackets in (42) then we obtain the expected maximum value. It should be noted that in general 

(max(.)) ^ max((.)). 

Proposition 2: Consider the instantaneous energy-efficiency of MIMO communication systems operating over ergodic time- 
varying i.i.d. Rayleigh channels with informed transmitter and receiver, and with instantaneous total receive power constraint. 
Then, the instantaneous maximum is achieved by the energy-efficient water-filling (EEWF) power allocation solution with 
per-channel variable power levels according to Proposition 1. The expected optimum energy-efficiency, the corresponding 
transmission rate and the total transmit power are given by 

rf = (rf , R a = (R x > , P t a = (P x , ) (44) 

where rf is the maximum of the instantaneous energy-efficiency and R x and P x are the corresponding instantaneous 
transmission rate and transmit power, respectively. The power allocation is performed dynamically, i.e., as per channel 
realization, according to a given total receive power constraint. 

This result is straightforward and follows from Proposition 1 above. 

In addition to (44) we also consider the average of the number of transmission channels 

N a = (N X ), (45) 

The transmission rate corresponding to the maximum energy-efficiency always satisfies R a < C e , where C e is the ergodic 
channel capacity of the time-varying i.i.d. Gaussian channel [11]. 

The power efficiency criteria with total receive power constraint for the SISO channel implies channel inversion. Hence, the 
condition for signal-to-noise ratio equality becomes P t = (j^)Pr (observe that in general ^ (j^))- For the time-varying 
SISO channel it can be shown that (^-) — > oo. However, the dynamic range of realistic channels is seldom as low as null 
and is often only a few tens or hundreds times larger than the sample average. Moreover, a better form of power control is 
obtained by only compensating for fading above a certain cutoff fade depth known as truncated channel inversion [12]. So we 
will assume that 1 < (-^-) < oo for our Rayleigh fading channels. 

Corollary 2.1: Consider a energy-efficient in the sense of Proposition 2 time-varying SISO channel. Then, the ratio between 
the transmission rate corresponding to the energy-efficient transmission R a and the ergodic SISO capacity C e satisfy the 
inequalities 

1 7? a 

The corresponding bounds for the ratio between the optimum SISO energy-efficiency rf and the energy-efficiency corre- 
sponding to the capacity achieving case rf, are 

rf 1 

1 - w - ( r>' (47) 

where we have assumed equal transmit signal-to-noise ratio such that P t = (^-)P r , 1 < (j--) < oo and (Ai) = 1. The upper 
and lower bounds in (46) and (47) are achieved at the high and low signal-to-noise ratio regimes, respectively. 
The derivations are given in Appendix B. 

It is worthwhile to note that if we have an ideal SISO Rayleigh channel then (^-) — > oo and the (average) total transmit 
power will therefore also increase infinitely P t a — > oo. 

IV. Numerical Examples 

In this section we present numerical computations that illustrate the dependence of the average of the optimal energy- 
efficiency on the signal-to-noise ratio and the number of antennas according to Proposition 2 above. We base our computations 
on 10000 realizations of the entries of the MIMO channel matrix H generated according to the assumption of i.i.d. Rayleigh 
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Fig. 1. (a) Average energy-efficiency v.s. signal-to-noise ratio; the different curves correspond to a different number of antennas N £ {1, 2, 4, 8, 16} going 
from the bottom to the top of the plot in increasing order and (b) average energy-efficiency v.s. number of antennas; the different curves correspond to different 
signal-to-noise ratios SNRf € {0, 5, 10, 15, 20} dB going from the bottom to the top in increasing order. The continuous lines correspond to EEWF and 
the dashed lines correspond to WF. 

(a) (b) 




SiWJJ, [clD] N 



Fig. 2. (a) Average transmission rate corresponding to Figure 1 v.s. signal-to-noise ratio; the different curves correspond to a different number of antennas 
N 6 {1, 2, 4, 8, 16} going from the bottom to the top of the plot in increasing order and (b) average energy-efficiency v.s. number of antennas; the different 
curves correspond to different signal-to-noise ratios SNR® £ {0, 5, 10, 15, 20} dB going from the bottom to the top in increasing order. The continuous 
lines correspond to EEWF and the dashed lines correspond to WF. 



(a) (b) 




SJViSJ, [dB] N 



Fig. 3. a) Average total transmit power corresponding to Figure 1 v.s. signal-to-noise ratio; the different curves correspond to a different number of antennas 
N 6 {1,2,4,8, 16} going from the top to the bottom the plot in increasing order and (b) average energy-efficiency v.s. number of antennas; the curves 
coincide for different signal-to-noise ratios SNR^ £ {0, 5, 10, 15, 20}. The continuous lines correspond to EEWF and the dashed lines correspond to WF. 



(a) (b) 




4 8 12 16 20 1 2 4 6 8 10 12 14 16 
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Fig. 4. (a) Average number of transmission channels corresponding to Figure 1 v.s. signal-to-noise ratio; the different curves correspond to a different number 
of antennas N £ {1, 2, 4, 8, 16} going from the bottom to the top of the plot in increasing order and (b) average energy-efficiency v.s. number of antennas; 
the different curves correspond to different signal-to-noise ratios SNR® £ {0, 5, 10, 15, 20} dB going from the bottom to the top in increasing order. The 
continuous lines correspond to EEWF and the dashed lines correspond to WF. 
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fading. Equal number of receive and transmit antennas TV has been assumed with channel normalization (\Hij\ 2 ) = 1. The 



Throughout the (a)-plots, from Figure 1 to Figure 4, we use the following labeling. The continuous lines represent the 
optimum energy-efficiency solution, i.e., the EEWF algorithm with total receive power constraint P r = 1, while the dashed 
lines correspond to the power allocation achieving the ergodic MIMO capacity, i.e., the WF algorithm with total receive power 
constraint P t = 1. Furthermore, in Figure 1(a), Figure 2(a) and Figure 4(a) the different curves correspond to a different number 
of antennas TV 6 {1,2,4,8, 16} going from the bottom to the top of the plot in increasing order. In Figure 3(a) the different 
curves also correspond to a different number of antennas but instead the count is going from the top to the bottom of the 
plot in increasing order. In Figure 1(b), Figure 2(b) and Figure 4(b) the different curves correspond to different signal-to-noise 
ratios SNR^ £ {0, 5, 10, 15, 20} dB also here going from the bottom to the top in increasing order. 

Figure 1(a) and Figure 1(b) show the energy-efficiency 77° as a function of SNR? and TV, respectively. As we can see, 
the energy-efficiency increases as a function of both STVi?" and TV. Increasing the number of antennas at a fixed SNRf 
results in a larger energy-efficiency as compared to increasing the signal-to-noise ratio at a fixed TV. Moreover, the energy- 
efficiency corresponding to the EEWF algorithm can be more than an order of magnitude greater than the energy-efficiency 
corresponding to the ergodic channel capacity. For example, in Figure 1(a), compare the top continuous line (EEWF) with 
the top dashed line (WF) for let's say SNR® = 20 dB; in both cases TV = 16. However, for the SISO case, the energy- 
efficiency corresponding to the WF power allocation is larger than the EEWF power allocation, which is contrary to the 
statement of Proposition 2 and more specifically of Corollary 2.1. Indeed, as we mentioned above, we have assumed that 
P r = 1 and P t = 1. Hence, they do not satisfy one of the conditions of Corollary 2.1 since (j-) ^ 1. In this case we have 
that if/Vc'e = lo g2(! + <STV.R t 7(-i-))/(log 2 (l + XiSNRf)) < 1. This follows from (63) as shown in Appendix B, where 
we have assumed P r = 1 for the energy-efficiency maximization and P t = 1 for the system capacity calculation. Hence, the 
signal-to-noise ratio SNR® is given by (xj-)/c>o an d V°o> respectively. 

For the MIMO case, i.e., TV > 1, the optimality of the EEWF power allocation comes at the expense of being suboptimal in 
terms of transmission rate. This is illustrated in Figure 2(a) and Figure 2(b) where the average transmission rate R a is plotted as 
a function of SNRf and TV, respectively. The dashed lines denote the ergodic MIMO capacity. As expected, the transmission 
rate corresponding to the EEWF algorithm is bounded from above by the channel capacity. The difference increases with 
both the number of antennas and the signal-to-noise ratio. The latter means that the total transmit power corresponding to the 
EEWF must be lesser than or equal to the total transmit power corresponding to the WF. This is shown in Figure 3(a) and 
Figure 3(b), where the total transmit power P t a is shown as a function of SNR® and TV, respectively. On one hand, it follows 
from Figure 3(a) that the total transmit power obtained does not depend on the SNR^. On the other hand, the total transmit 
power corresponding to the EEWF decreases with the number of antennas while the total receive power is fixed. Obviously, 
the total transmit power corresponding to WF shall not change with the number of antennas since it is a given constraint. 

Figure 4(a) and Figure 4(b) show the average number of transmission channels TV a as a function of SNRf and TV, 
respectively. As we can see, the number of transmission channels used by the EEWF approaches asymptotically the actual 
number of antennas as the signal-to-noise ratio increases. However, the convergence speed decreases with the size of the MIMO 
arrays. The WF algorithm uses all the available channels under the i.i.d. Rayleigh fading assumption. For the EEWF algorithm 
this is only valid for the SISO channel. 

As noted above, the maximum energy-efficiency provides a compromise between transmission rate and total transmit power. 
The EEWF power allocation results in reduced total transmit power at the same time as the achievable transmission rate 
experiments a loss with respect to the ergodic channel capacity. This trade-off should be taken into account when designing a 
wireless systems that aims at providing satisfactory performance in terms of both spectral- and energy-efficiency. 



The energy-efficiency of multiple-input-multiple-output (MIMO) systems with constrained received power establishes a 
useful trade-off between the transmission rate and the total transmit power. The proposed energy-efficient water-filling (EEWF) 
algorithm provides the optimal power allocation policy that maximizes the instantaneous ratio between transmission rate and 
the total transmit power given a minimum required total receive power in the single-user MIMO case. Studies of the static 
and the uncorrected fast fading Rayleigh channels show that the EEWF provides a non-zero optimal transmit power and 
non-zero optimal transmission rate. The optimal energy-efficiency as well as the corresponding transmission rate increase with 
the number of antennas and the signal-to-noise ratio. On the other hand, the total transmit power decreases with the number 
of antennas at the same time as the receive power requirement is satisfied. Future work may include channel correlation as 
well as the extension to multi-user MIMO scenarios. 



average signal-to-noise ratio is defined as SNRf 



-jr, where P t a is the average transmit power (49). 



V. Conclusions 



Appendix A 



We can, after some algebraic manipulations, arrive at the following modification of (22) 




(48) 
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where 

A, 

The variable x is related to a as follows 

A< = ln(2)Pf(na2 + P r )- (50) 

We can further recast (48) 

n. n n 

H(x- ai ) + J2 II (*-aj)=0- (51) 

»=i i=1 i=i J'/* 

Thus, solving (48) is equivalent to finding the roots of the monic polynomial of degree n under the assumption that x / ctj. 
Now we can expand the product terms into sums 

n n n— 1 

Pn (x) = y akxk + bMxk 

fc=0 i=l fc=0 

n—1/ n \ 



Y [ a k + Y hi ] x k + a n x n 

k=0 V i=l / 

Y ^ C 52 ) 



fc=0 



where 



c„ = 1, (53) 

n 

a- = Qfc + Y bki, k = 0, . . . n - 1. (54) 



The coefficients at are related to coefficients Ui through Vieta's formulas 



On-l 



O-n-2 



Y a ^ ^ 



i=l 



^«iOj, (56) 

l<i<j<n 

(-l) m ]T a 4l a l2 ■••«,„, (57) 

l<ii<i 2 <-<i m <n 



(-l)"II ai ' (58) 

i=l 

where the coefficient a„_ m is related to a signed sum of all possible sub-products of coefficients a»j m , taken m-at-a-time also 
known as elementary symmetric sums. We can apply the same relationships to obtain the coefficients bki from the coefficients 
on such that i ^ k. 

Invoking the Fundamental Theorem of Algebra we have that there are n complex roots that solve (50). However, we have to 
satisfy the condition that /i is a positive real number. Hence, we choose the solutions that are real positive x > 0. Furthermore, 
if there are several candidates we need to choose fi that actually maximizes the energy-efficiency. Let's assume that we have 
two solutions Hi and fi 2 with corresponding energy-efficiencies and transmit powers rjf, Pf x and rfe, P t 2, respectively. From 
the conditions of the problem we have that P r \ = P r2 = P r - Using this equality together with (22) and after discarding the 
noise term since it is asummed equal in both cases we obtain 

^ •(tf-AiiiSAOfof-^Ai) ' 

We can set the numerator to zero for each term of the sum and since > 0, — /iiP^A^) > and (tj| — ^Pn^i) > 
we obtain 

V2 -rff-Xi (M2^ 2 - Miifi) = 0. (60) 

Now, if r) 2 > r/f and P'f 2 < Pfi according to the premises of our problem we obtain that [i 2 > Hi must be satisfied. Hence, 
the largest positive real /i gives the solution to our problem if it exists. 
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Appendix B 

Let's first give the transmission rate to the capacity ratio 

R* l0g 2 (l + 



C e (lo g2 (l + ^))' 
and the corresponding energy-efficiency ratio 

log 2 (l ' 



(61) 



In order to show the right hand side of (46) and (47) we need to show that 



where 



P A P 

log 2 (l + 7T -^)<<log 2 (l + ^)), (63) 



(I)-i = — (64) 



is the geometric mean (GM) of the discrete realizations of A, i.e., Aj > 0, i = 1, . . . ,n. Let's consider the harmonic mean- 
geometric mean (HM-GM) inequality 



T n H _j_ <n(A i7 )", (65) 



where 7 = > is a constant. Hence 



iog 2 (i + ™ _jj < iog 2 (i + n^)" ). <&> 

n 

< iog 2 (n( i+A ^)")' ( 6? > 

since the logarithm function is a monotonically increasing function. The inequality in (67) follows from Mahler's inequality 
[13]. Finally, we readily arrive at 

log 2 (l + 7Tv) = log 2 (l+ 2 ) (68) 

n 

< log 2 n(!+^7)" (69) 

i=i 

1 " 

= - Vlog 2 (l + A l7 ) (70) 
n ^-^ 

i=l 

= (log 2 (l + A 7 )), (71) 

The equality is achieved for infinitely large signal-to-noise ratio, i.e., — > oo, which is obtained straightforwardly by applying 

CT ° p 

the L'Hospital's Rule to (46), where the derivative is taken with respect to 
To show the left hand side of (46) and (47) we need first to show that 

log 2 (l + 7#^)>7ivlog2(l + 4)- ( 72 ) 

But this is readily identified as the Bernoulli inequality, where y^r < 1. Now, we see that since the logarithmic function is 
concave function the Jensen's inequality for the expected value reads as 

(log 2 (l + ^))<log 2 (l + ^). (73) 

Hence, combining (72) and (73) into (61) and (62) gives the sought-after result. In this case, the equality is obtained for 
infinitely small signal-to-noise ratio by following a procedure similar the described above. 
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